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In the heavy quark effective theory, hadronic matrix elements of currents between 
two hadrons containing a heavy quark are expanded in inverse powers of the heavy 
quark masses, with coefficients that are functions of the kinematic variable v ■ v ' . For 
the ground state pseudoscalar and vector mesons, this expansion is constructed at 
order l/rug. A minimal set of universal form factors is defined in terms of matrix 
elements of higher dimension operators in the effective theory. The zero recoil nor- 
malization conditions following from vector current conservation are derived. Several 
phenomenological applications of the general results are discussed in detail. It is ar- 
gued that at zero recoil the semileptonic decay rates for B — > D I v and B — > D*£ v 
receive only small second order corrections, which are unlikely to exceed the level 
of a few percent. This supports the usefulness of the heavy quark expansion for a 
reliable determination of V c b- 
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I. INTRODUCTION 



Recent developments in the theory of heavy quarks have increased the prospects both 
for a reliable determination of some of the fundamental parameters of the standard model, 
and for a study of nonperturbative QCD in the weak decays of heavy mesons and baryons. 
The excitement is driven by the discovery of a spin-flavor symmetry for heavy quarks that 
QCD reveals in the limit where the quark mass uiq — > oo, in which certain properties of a 
hadron containing the heavy quark become independent of its mass and spin These 
symmetries are responsible for restrictive relations among weak decay amplitudes and reduce 
the number of independent form factors. The description of semileptonic transitions between 
two ground state heavy mesons 0,|3J or baryons |3],^,||[7] becomes particularly simple. In the 
limit where the heavy quark masses are much larger than any other hadronic scale in the 
process, the large set of hadronic form factors is reduced to a single universal function of 
the kinematic variable v ■ v', which measures the change of velocities that the heavy hadrons 
undergo during the transition. It depends on the quantum numbers of the light degrees of 
freedom, but not on the heavy quark masses and spins. In addition, the conservation of 
the vector current implies that this celebrated Isgur-Wise form factor is normalized at zero 
recoil, allowing model-independent predictions unaffected by hadronic uncertainties. 

Clearly, a careful analysis of at least the leading symmetry-breaking corrections is es- 
sential for any phenomenological application of the heavy quark symmetries. An elegant 
framework in which to analyze such corrections is provided by the so-called heavy quark 
effective theory, which allows for a systematic expansion of decay amplitudes in powers of 
1/rn.Q [p|,p|, p!0| , pT| , p!2| , p!3| , p!4| . The coefficients in this expansion are given by matrix elements 
of operators in the effective theory and can be parameterized in terms of universal form 
factors, which characterize the properties of the light degrees of freedom in the background 
of the static color source provided by the heavy quark. At leading order one recovers the 
Isgur-Wise limit, in which only a single function remains. But already at order I/trq one 
encounters a larger set of universal form factors, which affect all but very few of the sym- 



metry predictions that hold in the infinite quark mass limit |IqJlqJ17|l . An understanding 
of these functions is at the heart of nonperturbative QCD, but it is ultimately necessary 
for any quantitative analysis based on heavy quark symmetries. In the future, one might 
hope to compute the universal form factors from first principles by using a formulation of 
the effective theory on a lattice ||,pl|,|l9[| . In the meantime, QCD sum rules ||2(J offer a less 
ambitious approach to this problem, and have recently been employed to study the decay 
constants of heavy mesons, the Isgur-Wise form factor, and the universal functions that 
appear at order 1/toq in the heavy quark expansion [ 17 |H , 23 HE \ . One can also gain 
valuable information about symmetry-breaking corrections from measurements of certain 
ratios of form factors | |25| | . 

In this paper we analyze current- induced transitions between ground state heavy mesons 
at order 1/ttiq in the heavy quark expansion. Such an analysis is particularly relevant 
for the important cases where the leading 1/rriQ corrections are known to vanish at zero 
recoil. An example is the semileptonic decay B — > D*£u, which therefore seems ideal 
for a measurement of the weak mixing parameter V c b PiP6|- In Sec. |TI| we discuss the 
parameters of the effective theory that appear at subleading order. The general method of 
computing power corrections is outlined in Sec. [TTT| , together with a review of the analysis 
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of the 1/toq corrections to transitions between heavy mesons. In Sec. |IV| we extend this 
analysis to second order. We identify a minimal set of universal functions and give their 
relation to matrix elements of higher dimension operators in the effective theory. The zero 
recoil normalization conditions imposed on some of these form factors are derived. Although 
in principle straightforward, the analysis is tedious and involves considerable technicalities of 
the heavy quark effective theory. The reader not interested in these details is encouraged to 
proceed to Sec. |V|, where we summarize our results and illustrate them for some specific cases 
of phenomenological relevance. In particular, the corrections affecting the determination of 
V c b from exclusive semileptonic B decays are investigated. We also study the fictitious limit 
of vanishing chromo-magnetic interaction, which leads to great simplifications and might 
serve as an estimate of the dominant corrections. 

Based on the analysis for heavy mesons, the 1/wq corrections to matrix elements between 
heavy baryons can readily be derived. We discuss this subject in the following paper p7| . 



II. PARAMETERS OF THE EFFECTIVE THEORY 

The construction of the heavy quark effective theory (HQET) is based on the observation 
that, in the limit tuq ^> Aqcd, the velocity of a heavy quark is conserved with respect to soft 
processes. It is then possible to remove the mass- dependent piece of the momentum operator 
by a field redefinition. To this end, one introduces a field hq(v , x), which annihilates a heavy 
quark with velocity v a , by 



h Q (v,x) = e im ^ x P + (v)Q(x), (2.1) 

where P+{v) = |(1+^) is an on-shell projection operator, and Q(x) denotes the conventional 
quark field in QCD. If P a is the total momentum of the heavy quark, the new field carries 
the residual momentum k a = P a — m,QV a . 

There is obviously some ambiguity associated with the construction of HQET, since the 
heavy quark mass used in the definition of the field hq is not uniquely defined. In fact, for 
HQET to be consistent it is only necessary that k a be of order of Aqcd, stay finite in 
the limit itlq — > oo. It is intuitively clear that different choices for rriQ must lead to the 
same answer for any physical matrix element, and this can indeed be shown to be the case 
2"8fl . Yet it is advantageous to adopt a special choice for which the resulting effective theory 



becomes particularly simple, in the sense that there are no "residual mass terms" for the 
heavy quark and the heavy quark expansion becomes a covariant derivative expansion. This 
prescription provides a nonperturbative definition of the heavy quark mass, which has been 
adopted implicitly in most previous analyses based on HQET. It is important to realize, 
however, that so defined, the mass uiq is a nontrivial parameter of the effective theory. 
In the limit vhq — > oo, the effective Lagrangian for the strong interactions of the heavy 



quark becomes p"T|,|T2"|.p"3 



£hqet = h Q iv-DfiQ , (2.2) 

where D a = d° — ig s t a A^ is the gauge-covariant derivative. For finite tuq, there appears in 
the Lagrangian an infinite series of power corrections involving higher dimension operators, 
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£power — ~ £l + ~ T £-2 + ■ ■ ■ ■ (2-3) 

2mn Am 2 n 



Note that it is natural to expand in powers of l/2mg since, after the field redefinition ( |2.1|) , 
2rriQ is the mass associated with the heavy antiquark field which is integrated out ||13|| . 
Omitting an operator whose matrix elements vanish by the equation of motion, the leading 
term in (|2.3| ) is given by |14| 



Ci = h Q {iDfhq + Z(m Q /fi) h Q s aP G a/3 h Q , (2.4) 

where s a p = —\<J a @-, and G a/3 = [iD a , iD 13 } = igstaG^ 13 is the gluon field strength. In leading 
logarithmic approximation, the renormalization factor for the chromo-magnetic moment 
operator is 



Z{m Q /fi) 



H9//3 



; P = 33-2n f , (2.5) 



where rif is the number of light quark flavors with mass below tuq. The kinetic term in 
( f2.4| ) is not renormalized. 



The purpose of the heavy quark expansion is to make the mQ-dependence of some 
hadronic quantity A explicit by writing 

A{m Q ) = C (m Q /fi) A (fi) + -— C x {m Q /^) A 1 ( f i) + • ■ ■ , (2.6) 

in such a way that the coefficients Ai{fi) are universal, mQ-independent parameters, and 
C%{v(IqI ' \x) are purely perturbative coefficients, which dependent on rng only via the running 
of the strong coupling a s (mg). The aim is to relate Ai to matrix elements of operators 
in HQET evaluated between the eigenstates of the lowest order Lagrangian £hqet- This 
paper focuses on the ground state pseudoscalar and vector mesons, which form a degenerate 
doublet under the heavy quark spin symmetry. These mesons have the same velocity as the 
heavy quark which they contain. Their common mass M, however, differs from the mass 
of the heavy quarks by a finite amount A = M — tuq, which measures the "mass" carried 
by the light degrees of freedom. Because of the field redefinition ( p.l|) , it is this mass which 
governs the x-dependence of states in the effective theory: 

|M(x)) HQ et = e- iA ^|M(0)) HQ ET • (2.7) 

A is a universal parameter which can be defined in terms of a matrix element of a higher 
dimension operator in HQET. Using the equation of motion iv-D Hq = 0, which follows from 



the effective Lagrangian £hqet, it is easy to see that |28 



(0\gtvDTh Q \M(v)) 

(0\qTh Q \M(v)} • ( - } 

Here T is an appropriate Dirac matrix such that the currents interpolate the heavy meson 
M. This relation shows that A is in fact a parameter describing the properties of the light 
degrees of freedom in the background of the static color source provided by the heavy quark. 
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It turns out that this mass scale also enters the leading power corrections to heavy meson 
form factors and determines the canonical size of deviations from the infinite quark mass 



limit |15| ) |16| . A recent analysis of A using QCD sum rules predicts [|i7f 

A = 0.50 ± 0.07 GeV. (2.9) 

The eigenstates of £hqet differ from the states of the full theory. In particular, their 
mass M differs from the physical masses of pseudoscalar or vector mesons by an amount of 
order 1/mq. These mass shifts are computable in HQET. The physical masses itlm obey 
a heavy quark expansion, which we write as (rriM — wq) = A + Am 2 M /2mQ + - • - . In the 
meson rest frame, 

(M W |( A)|M W ) 
(M(v)\hyi Q \M(v)) 

A convenient way to evaluate hadronic matrix elements in HQET is by associating spin wave 
functions 

M(v) = VMP + (v) ( T 75 ; P se + udoscalar ™ P (2.11) 
[f ; vector meson V 

with the eigenstates of £hqet P, [29| , |30| . These wave functions have the correct transforma- 
tion properties under boosts and heavy quark spin rotations. Here e a denotes the polariza- 
tion vector of the vector meson. For reasons of simplicity we shall often omit the argument v 
in both P + and A4. We note that Ai = P + M. P_, where P± = |(1 ± f). Lorentz invariance 
allows one to write any matrix element as a trace over these wave functions and appropriate 
Dirac matrices. For the matrix elements in ( |2.10| ) we define hadronic parameters Aj by 



(M\ h Q {iDfh Q \M) = -Ai tr{ M M } = 2M Ai , 
(M| hqs^G^hq \M) = -X 2 (ji) tr{ icr af} Ms a/3 M } = 2d M M X 2 (fi) , (2.12) 

where dp = 3 for a pseudoscalar meson, and dy — — 1 for a vector meson. The conservation 
of the vector current implies that, in the rest frame, the matrix element in the denominator 
is given by (M\ hqhq \M) = 2M. We thus have 

Am 2 M = -Ax - d M Z(n) A 2 (/i) . (2.13) 

The universal parameters Ai and A2 are the analogs of A at subleading order in the heavy 
quark expansion. They are independent of mq. Whereas Ai is not renormalized, ^(l^) 
depends on the renormalization scale in such a way that the product Z(fi) A2(jw) is scale- 
independent. 

An estimate of the value of A2 can be obtained from the measured mass splitting between 
the B* and B mesons, assuming that higher order corrections in the B system are small. 
One finds 

m|. - m| w Am|, - Am| = 4A 2 (m fe ) w 0.48 GeV 2 , (2.14) 



where the experimental value has been taken from Refs. [31]. Using (|2.5| ) for the evolution 



of this parameter down to the low energy scale 2A w 1 GeV, we obtain 
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A 2 (2A) « 0.15 GeV 2 



(2.15) 



Unfortunately, it is not possible directly to relate the spin-symmetry conserving parameter 
Xi to an observable quantity. Recently, QCD sum rules have been used to compute both 
Ai and A 2 [I7|. The spin-symmetry breaking correction was found in excellent agreement 
with experiment, A| n = 0.12 ± 0.02 GeV 2 , and a rather large value for the spin-symmetry 
conserving correction was obtained, Xf T - ~ 1 GeV 2 . However, the sum rule analysis suggests 
that it might be more appropriate to use an effective value of Ai in the b and c system which 
could be substantially smaller, even compatible with zero. A measurement of Ai on a lattice 
could help to clarify this issue. 



III. MESON FORM FACTORS IN THE EFFECTIVE THEORY 

Let us now review the analysis of current-induced transitions between two heavy mesons 



to subleading order in HQET, as performed by Luke |15| . This will help to outline the general 
procedure and set up the conventions we will need in Sec. [TV]. The aim is to construct the 
heavy quark expansion ( |2.6| ) for matrix elements of the type (M'(v')\ Q'TQ \M(v)), where 
T is an arbitrary Dirac matrix. In this case the universal parameters are functions of the 
kinematic variable w = v ■ v ', and the perturbative coefficients, subsequently denoted by 
Cj,Cj and c'j, depend on w and both heavy quark masses. The current Q'TQ has a short 
distance expansion in terms of operators of the effective theory. It reads 

Q'TQ ^J2 C jh' Fjh 

j 

+ ^Y, c jh' T P D ^ h + 7^—Y, c ' j h'(-'Da)T' a h + --- } (3.1) 

2m Q j J 2m Q > j J J 

where we have abbreviated h = hq(v) and h' = hqi{v'). The matrices Tj are in general 
different from V and can depend on v and v'. At tree level, however, one has 



E^Tj-r, ^ Cj . r ^r 7 Q , y>;r; a ^ 7 Q r. (3.2) 

j j j 

Using the trace formalism described in Sec. |TJ, matrix elements of the leading term in 



( |3.1D can be parameterized as 

(M'\h'Th\M) = -£(w,fi)tr{M'rM}, (3.3) 

where we omit the velocity argument in the states and the wave functions in order to simplify 
notation. It is to be understood that quantities without a prime refer to the initial state 
meson M, while primed quantities refer to the final state meson M'. Also, from now on 
m will designate a generic heavy quark mass. In general the form factor £ could be some 
matrix-valued function of v and v', but in this case the projection operators contained in the 



spin wave functions restrict it to a scalar function of w. Eq. ( |3.3|) implies that, to leading 
order in the heavy quark expansion, all matrix elements of currents between pseudoscalar or 
vector mesons are described by a single form factor, the Isgur-Wise function [|2|,|3|j29|] . The 
kinematical information is contained in the trace over spin wave functions. By evaluating 
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the special case of mesons with equal mass and velocity, one readily derives the zero recoil 
normalization condition = 1 as a consequence of the conservation of the vector 

current. 

At subleading order in fl3.1| ) one encounters current operators which contain a covariant 
derivative. Their matrix elements are represented by the diagrams shown in Fig. |l|(a) and 
can be parameterized as 

(M'\h'T a tD a h\M) = -ti{^ a (v,v' } fi)M'T a M}, 
{M'\h! (-il5 a )Y' a h\M) = -tr{t a {v',v,v)M'r' a M}. (3.4) 

Note the interchange of the velocities in the second matrix element. The most general 
decomposition of the universal form factor £ Q involves three scalar functions. Following 
Ref. JL5J, we define 

f a (tf, V', fl) = £+(w, fl) (V + V') a + fl) (V - V') a - f 3 (ll7, fl) 7« . (3.5) 

T-invariance of the strong interactions requires that these scalar functions be real. Using 
( p.7|) and the fact that id a (h! Y h) =h'i D a T h + li'Y iD a h, one finds that 

£-(w,fi) = ^£(w,fi) . (3.6) 

This is where the parameter A enters the analysis. 

The equation of motion, iv ■ Dh = 0, yields an additional relation among the scalar 
form factors. Taking into account that under the trace £ Q is sandwiched between projection 
operators, one obtains 

P-v a t a (v,v',iJL)PL = 0. (3.7) 

For the remainder of this paper we use the symbol "=" for relations such as this, which 
are true when sandwiched between the projection operators provided by the meson wave 
functions. We thus write v a £, a {v , v', fx) = 0. In terms of the scalar functions this is equivalent 
to 

(w + 1) y)-(w- 1) fj) + U*>, V) = • (3-8) 

We shall use this equation to eliminate In particular, it follows that at zero recoil 
2£ + (l,/i) +£3(1,//) = 0. This relation has an interesting consequence, since it implies that 

Uv, v, = [2 £+(1, /i) + 6(1, /x)] v a = , (3.9) 



showing that matrix elements of the higher dimension currents in (3.1) vanish at zero recoil. 
This is the first part of Luke's theorem [TjJ. In its above form it is obvious that this result 
is true to all orders in perturbation theory p2| , since it does not rely on the structure of the 
perturbative coefficients in ( |3.1|) . 

A second class of 1/m corrections comes from the presence of higher dimension operators 
in the effective Lagrangian. Insertions of operators of C\ in ( p.3|) into matrix elements of the 
leading order currents represent corrections to the wave functions, which appear since the 



7 



eigenstates of £hqet are different from the eigenstates of the full theory. The corresponding 
diagrams are shown in Fig. 0(b). The relevant matrix elements can be written as 



(M'\i dxT{ J(0),£i(x)}|Af) 



(M'\i dxT{ J(0),£[(x)}\M) 



-A x {w,n) ti{ M'TM } 

-Z(m Q /fi) tr{ A aP (v, v', fi) M' T P + s a3 M } 

-Ax(w,n) tr{M'TM } 

-Z(mcy/(i) tr{ A al3 (v', v, //) M' s aP P' + V M } . 



(3.10) 



where J = h'Th is a lowest order current. Noting that v a P + s af3 Ai = 0, we write the 
decompositioni 



A a p(v, v', /i) = A 2 (w, fj) (v' a 7p - v'p'ja) + A 3 (w, fx) ia af3 . 



(3.11) 



The four independent functions £3 and Ai, as well as the mass parameter A, suffice to 
describe the first order power corrections to any matrix element of a heavy quark current 
between ground state mesons. To get a picture of the structure of the corrections let us for 
simplicity neglect radiative corrections. In this case, there is a simple relation between the 
currents in HQET and the current in the full theory. Consider now the power corrections 
proportional to l/niq. They leave the wave function of the final state meson unaffected, but 
change the simple structure of M.(v ). The part proportional to the on-shell projection oper- 
ator P + will be modified, and a component proportional to P_ will be induced, representing 
the "small component" of the full wave function. Hence 



M(v) -> P+(v) L^{v,v') + P_(v) L M (v,v') . 
The general form of L± is 

L^{v i v') = VM(- l5 )L 1 {w) } 



(3.12) 



■v 1 



L + \v,v\ 
L P (v,v>) 
L v _(v,v') 



<f.L 2 {w) + e-v' L 3 (w) 



y/M (-7b)L 4 H, 

fL B {w) + e-v'L G (w) 



(3.13) 



The insertions of higher order terms from the effective Lagrangian in ( |3.1U| ) obviously con- 
tribute to only. On the other hand, in the absence of radiative corrections the matrix 
elements of the higher dimension currents can be written in the form 

(M'\h'Tiph\M) = -tr{M' V P_[ la M C(v,v')}} , (3.14) 

where we have used (|3.7|) to insert P_ between T and j a . Consequently, these corrections 
contribute to L^ 1 only. This is expected since ip h is proportional to the small component 
of the full heavy quark spinor. By evaluating the relevant traces one easily obtains 



1 Our functions are related to those defined in Ref. [15| by Ai = 2\i, A 2 = —2x2 and ^3 = 4%, 
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L x = Ai + 2{w - 1)A 2 + 3A 3 



L 2 
Li 



Ax-A 3 , 
-2A 2 , 

-AC, 

2 (Ae+6) 



(3.15) 



w + 1 

and the complete matrix element becomes 

(M'\Q'TQ\M) = -t(w)tr{M'TM} 

— trf M'V \P + L A J(v,v')+P-L M (v,v')}} 

2m Q ^ JJ 



(3.16) 



2m Q/ 



— tr{ [l5 V, u) P; + L M '(w', v)P'_]tm} + 



This way of organizing the corrections reduces to a minimum the effort required to compute 
the traces. If radiative corrections are taken into account, it still suffices to define these 
six functions Lj, as long as one stays with the leading logarithmic approximation for the 



perturbative coefficients. The corresponding expressions are given in Ref. | 25fl , where the 
functions Lj were called ft. 

Let us evaluate fl3.16j ) for the matrix elements of the vector and axial vector currents, 
and An, between bottom and charm mesons, which can be described completely in terms of 
fourteen meson form factors hAw). We define 



(D{v')\ Vn \B(v)) = y^^B h+(w) (v + v% + h-(w) (v - v'\ 



(D*{v', e') | Vn \B{v)) = iy/m B m D * h v {w) e^p v' a v p , 



(D*{v', e 1 ) \A tl \B(v)) = y/m B m D . h Al (w){w + l)e 



-h A2 (w) e'*Wn - h As (w) e'*vv' M 



(3.17) 



{D*{v\ e') \ Vfj, \B*{v, e)) = y/m B *m D * 



e-e 



h x {w) (v + v')^ + h 2 (w) (v - v') f 



+h 3 (w) e'*-ven + h 4 (w)e-v'e'; 



i /* 
—e-v e -v 



h 5 {w)v tl + he(w)v' 



(D*(v', e') | An \B*(v, e)) = iy/m B .m D . e a e'* 13 [h 7 (w) (v + v') v + h 8 (v- v') 1 
At leading order in the heavy quark expansion one finds that 

h + = hy = fiAi = h As — hi — h% — /i4 = hi — £ , (3.18) 
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while the remaining six form factors vanish. The expressions arising at subleading order are 
given in Appendix B. Here we restrict ourselves to three important cases, namely 



(3.19) 



2mh 



w + 1 



+ 



2m r 



L 2 {w) -^—^-L 5 (w) 
w + 1 



The conservation of the vector current in the limit m& = m c implies the zero recoil normal- 
ization conditions h + (l) = hi(l) = 1, from which it follows that £i(l) = L 2 (l) = 0, i.e. 



A 1 (l, f t)=A 3 (l, f i) = 



[3.20) 



This is the second part of Luke's theorem, which is again true to all orders in perturbation 
theory. It follows that A a p(v,v, fi) =0, so that the matrix elements in ( p,10|) vanish at zero 
recoil. 

In summary, Luke's theorem implies that the matrix elements which describe the first 
order power corrections in HQET vanish at zero recoil. It is important to realize that this 
does not imply that the meson form factors are unaffected by 1/m corrections p3| |. In 
fact, the theorem only applies for form factors which are not kinematically suppressed as 
v' — > v. Besides h + and h\ those are Ha 1 and h 7 . The form factor h^, which according to 
( gig ) is indeed seen to be unaffected by first order power corrections at zero recoil, plays 
an important role in the determination of from semileptonic decays ||26|| . It is one of the 
purposes of the next section to investigate the second order corrections to this form factor. 



IV. SECOND ORDER POWER CORRECTIONS 

The analysis of higher order corrections in HQET makes use of the same techniques 
as those developed above. The second order corrections can be divided into three classes: 
corrections to the current, corrections to the effective Lagrangian, and mixed corrections. 
We shall discuss each of them separately below. In order to keep the presentation as simple 
as possible we will often ignore radiative corrections; however, we will always make clear 
how they could be incorporated into our analysis. 

A. Second Order Corrections to the Current 

At tree level, the expansion of the heavy quark current reads [cf. (|3.1| )1 

Q'TQ -> h'Th + ^— h'Tiph + h'(-iB)Th 

2m Q 2m Q i 
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Am Q m QI 



h'(-ip)Tiph + --- . (4.1) 



On dimensional grounds, the operators appearing at second order are bilinear in the covariant 
derivative (recall that G afi = [iD a ,iD (3 ]). This remains true in the presence of radiative 
corrections, although a number of additional operators are induced. It thus suffices to 
consider the single hadronic matrix element 

(M'\h! {-iD a )T ap iD p h\M) = -ti{^{v,v',fi)M'T a/3 M}, (4.2) 

represented by the third diagram in Fig. 0(a). From now on we shall always omit the 
//-dependence of the universal form factors except in the equations which define them. Con- 
sidering the complex conjugate of the above matrix element, one finds that the form factor 
must obey the symmetry relation 

^PaW, V) = 1pap(v, v') , (4.3) 

which reduces the number of invariant functions to seven. It is convenient to perform a 
decomposition into symmetric and antisymmetric parts, i/i a p = ~ [i/j^p + V'asl > an d to define 

1P^(V, V') = W 9a(3 + ^O) {V + V') a (v + V')p + lP§(w) (V - V%{v - v')p 
+^i{w) [(« + v') a Jp + (V + ?/)/37oJ , 

(4.4) 

V'i(U) V> ) = ^ti W ) {VaV'p - v' a Vp) + l/4(w) [{V - f %7/3 - (V - v')p^ a 

As in ( |3.7| ) one can use the equation of motion to derive relations among the scalar form 
factors. It follows that under the trace 

v^ a p(v 7 v') = 0, v' a i; a p(v, v')=0. (4.5) 

These conditions are equivalent because of (|4.3|) and lead to the three relations 



+ (w + 1) ^ -(w-i) $ - i)f + wip? - ^ - ipf = o , 
(w + 1) + (w - 1) - i>! - + ^2 = o , 

(u; + l)Vf + (w-l)^-^ = 0, (4.6) 

which reduce the number of independent functions to four. 

One can use an integration by parts to relate (f4.2| ) to matrix elements of operators 
containing two derivatives acting on the same heavy quark field, which are represented by 
the first two diagrams in Fig. 0(a). It follows that 

(M'| h! T af3 iDJDp h\M) = -tr{ iff^v, v') Jd' T af3 M } 

—A (v - v') a tr{ ^(v, v') M' T a ?M } (4.7) 



with £p as defined in (|3.4|) . Matrix elements of operators with both derivatives acting to the 
left can be obtained in a similar way. In particular, we may derive from (|4.7| ) the matrix 
elements 
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(M'\ ti T (iDfh \M) = -0oH ti{M'TM}, 
(M'\ h! T af3 G a p h\M) = -tr{ 4> a p(v, v') M' T a(3 M } 



(4.8) 



the second of which is needed in 
we find 



If). Choosing the same decomposition for Q/ 3 as for i/r^, 



= 2i/>? + (w + l) ^ -(w-1) ^ - 2^f - A V - 1)£ . 

1 = t# + ^[A>-lK-2A£ 8 ], 



w+l 

3 



i2 = ^-AC: 



(4.9) 



We have already encountered matrix elements similar to ( |4.8|) in the discussion of mass shifts 
in Sec. ||, and from a comparison with (|2.12|) we find the zero recoil conditions 



^o(l) = A 1 , 
k(l) = A 2 , 



o Ax + 9 A^ 



(4.10) 



the last one being a consequence of the relations ( f4.6|) , which allow us also to express the 
form factors ipf in terms of the functions 0j. After some algebra we find 



Wl = 00 + W 01 

1 



w 



(202 + 3 ) 



w — 1 



w 



2(w + i; 
1 



2(w + iy 



w + l 
6 + (2w - 1) 0x] 

20 2 + (2w + 3) 03 - (2 - iu)(u; - 1) A 2 £ - 4(iu - 1) A & 



(4.11) 



+ 



1 



20 2 + 03 + 2wA 2 £ 



5 



W + l 

We have introduced the function 

1 



4(w + l) 

(w-l)0 2 + 03 -(™-l)A63 



0(w) 



It' 



o (w) + (w + 2) 1 (w) - 30 2 H-| 3 (w) 



(4.12) 



which is nonsingular as w — > 1 because of ( |4.10| ). 

The above relations allow us to prove a theorem which is the analog of the first part of 
Luke's theorem: 

Theorem 1: At zero recoil, matrix elements of second order currents 
in the heavy quark expansion can be expressed in terms of Ai and A 2 . 

For the proof we note that, to all orders in perturbation theory, the relevant operators contain 
two covariant derivatives. Because of ( |4.2| ) and ( |4.7| ) the corresponding matrix elements at 
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zero recoil only involve ip a p(v,v) sandwiched between projection operators. Using ( |4.10|) we 
find that 

il>ap(y, v) = [g a p - v a vp] y + ia a/3 y , (4.13) 

which proves the theorem. Furthermore, we note that at tree level only the last term in 
( O ) contributes at zero recoil, since vP<f> a p(v,v) =0. The corresponding corrections are of 
order \i/m Q m Q/ . 



B. Second Order Corrections to the Lagrangian 



Apart from operators whose matrix elements vanish by the equation of motion, the most 
general form of the coefficient L 2 appearing at second order in the expansion of the effective 
Lagrangian in (|2.3|) contains two terms, 



C 2 = Z 1 {m Q /n)hvpiD a G a/3 h + 2Z 2 {m Q / f i)hs a pv 1 tD a G (3 "<h. (4.14) 
In leading logarithmic approximation the renormalization factors are given by [ j34[ 



Zi(m Q fn) 



55 



Z 2 (m Q /n) = — - — 




(4.15) 



These operators have the same Dirac structure as the operators in C\ in fl2.4|) , and conse- 
quently their matrix elements are of the same form as those of C±. In analogy to ( |3.10| ) we 
thus define 

(M'\i fdxT{J(0),C 2 (x)}\M) : 



-Z x B^w, fi) tr{ M T M } (4.16) 
-Z 2 tr{ B aP {v, v', n) M' T P + s a/3 M } , 

and similarly for an insertion of C 2 . As before, J = h'Th denotes a lowest order current. 
The corresponding diagrams are the first two shown in Fig. 0(b) . The decomposition of B a p 
is of the same form as that for A a p in ( |3.11|) . It involves two functions B 2 and B 3 . 

Another type of 1/ m 2 corrections comes from a double insertion of the first order correc- 
tion C±, as shown in the third and fourth diagrams in Fig. 0(b). The corresponding matrix 
elements have a more complicated structure. We define 

(M'| % J dxdy T { J(0), £i(a;),£i(y) } |M> 

= -d(w, n) tr{ M' T M } - Z tr{ C a/3 (v, v', fj) M' T P+ s af3 M } 



-Z 2 tr{ C a(3l& (v, v', n) M' T P+ s a ?P + s~< s M } . 



(4.17) 



Again, the corresponding matrix elements with two insertions of £[ can be obtained by 
conjugating the matrix elements as in fl3.10| ). The decomposition of C a p is the same as 
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that for A a p, involving two form factors C 2 and C3. The most general decomposition of the 
four-index object C a ^ involves nine invariant functions, C 4 to C\2- They can be defined 
by 

C a pys(v, v') = C 4 (w) (g ai gps - gasgp-y) + C 5 (w) a lS a a p 

+ C 6 (w) {g al iep5 ~ gp^aS ~ 9aS^^ + 9 05™ art) 

+ C 7 (w) (v'^js - v' s 'yy)(v' a 'yp - v' p j a ) 

+ C 8 (w) (gaWpv' d - gpyV^v'g - g^v^ + gp&v'^) 

+ Cg(w) (g al v'f3l8 ~ g^v'alS ~ QagV^ + ^<7 7 ) 
+ C W (w) (gorylpv's ~ g^laV's ~ gaSlpV^ + flf/MTa^) 

+ C n (w) (ia a ~ / v'p'ys - ia^v'^g - ia aS v'^ + iap S v' a ^) 

+ C 12 (w) {ia a ^pv' s - iopfiav'i - ivasjpv'^ + iopslav'J ■ (4.18) 

Finally, there are corrections resulting from insertions of both C\ and £[, as shown in 
the last diagram in Fig. 0(b). They have the form 



(M'| i 2 j dxdy T { J(0), d(x), C'M } \M) 



= -D 1 (w,n)ti{M TM} 

-| tr{ D a/3 (v, v', y) M' T P + s af3 M } 

-y tr{ D aP (v', v, fj) M' s al3 P[ T M } 

-ZZ' tr{ D a ^ s {v, v', /i) M' s q/3 P; V P + s~< 5 M } . (4.19) 

The form factor D a p is again of the same form as A a p and involves two functions, D 2 and 
D3. The most general decomposition of the four-index object D a ^s is similar to that of 
C a /3 y s- However, because of the symmetry of the matrix element ( |4.19|) this quantity has to 
obey the constraint 

D a p y6 (v,v r ) = D ySa p(v',v), (4.20) 
which allows only seven independent functions, D4 to Diq. We choose the decomposition 

D a p j5 (v, v') = D 4 (w) (g a7 gps - 9as9i3 7 ) + D 5 (w) 0^0 - a(3 

+ D 7 (w) - vfaXVaTfp - V^ a ) 

+ D 8 (w) {g ai v p v' 5 - g/3jv a v' g - g^vpv'^ + gp S v a v' 7 ) 
+ D Q (w) g ai vpl6 - gp~fV a ~i s - gasvpj^ + gpsValj 

+ D 10 (w) v^bio a ~ t - Vcfl&iopn ~ v^ia aS + v^iaps 

+ i(? al lpv' & - icTfrJaV's - iVaSjpV 1 + IGps^v' . ( 4 -21) 
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In total, twenty-five universal functions Cj and Di are necessary to parameterize the 
effects of second order corrections to the effective Lagrangian of HQET. Unlike the correc- 
tions to the current, there are no relations imposed on these form factors by the equation of 
motion. 

Before proceeding, we have to discuss an additional source of second order corrections, 
which is related to the ones encountered above. As discussed in Sec. O, the mass M in the 
wave functions that we associate with the eigenstates of £hqet is different from the physical 
mass rriM- It is the physical mass, however, that appears in the normalization of matrix 
elements of the vector current, which one uses to derive zero recoil conditions for some of 
the universal form factors. At second order in the heavy quark expansion one has to take 
into account this difference and perform a mass renormalization of the wave function, 



M(v) 



z l/2 



J M 



M(v), 



7 l/2 
J M 



M 



in the first term in (|3.16|) . This is compensated by a counterterm 



7V2 7V2 



£{w)tT{M TM} 



'H + H:) £Wtr{ ^ rjw}) 



4mi 



Ami 



(4.22) 



(4.23) 



which, according to ( |2.13|) , effectively adds Ai £ to C\ and \ 2 i to C3. 



C. Combined Corrections to the Current and the Lagrangian 

The third and last type of 1/m 2 corrections arises from the combination of first order 
corrections both to the current and to the Lagrangian, as shown in Fig. 0(c). The relevant 
structures are 

(M'\iJ&xT{h'T 1 iD 1 h,C 1 {x) } \M) = -ti{E^(v,v',fj.) M't^M } 

-Z tr{ E^ p (v, v 1 , fi) M' r<P + s af3 M } , 

(4.24) 

(M'\ijdxT{ h! (-W 7 )Th,Ci(x) } \M) = -tr{ E'^v, v', fi) M'V^M } 

-Z tr{ E'^(v, v', fi) M' r<P + s a(3 M } . 

As previously, insertions of £[ give rise to the conjugate matrix elements. The form factor 
E 7 may be parameterized as 

E 7 (v, v') = Et(w) + E 2 (w) v' n + E 3 (w) 7 7 . (4.25) 

The most general decomposition of E ia p involves eight functions, which we define by 

E ia p(v, v') = (v' a jp - v'pja) E 4 (w) v 1 + E 5 (w) t> 7 + E 6 (w) 7 7 



+icr a p E 7 (w) v 1 + E 8 (w) v'^ + Eg(w) 7 7 
+{g ai E w (w) v'p + E n (w) 7/3 - (a <-+ 
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0)} 



(4.26) 



The equation of motion implies v 7 -E 7 = and v 1 E ia a = 0, which is equivalent to 



E 1 + wE 2 - E 3 = 0, 
£ 4 + w E 5 + E 6 = , 
E 7 + w E 8 - E 9 = . 



(4.27) 



Here we have used the fact that v a P + s a @A4 = 0. 

We define form factors E[{w) by identical decompositions. In this case, the equation of 
motion leads to the relations 



wE[ + E , 2 -E' 3 = 0, 
wE', + E' 5 -E' 6 + E' n = 0, 
w E' 7 + E' 8 - E' g = . 



(4.28) 



These functions are not independent of Ei, however, the reason being that the matrix ele- 
ments in ( 4.24| ) are related to each other by an integration by parts. This relation has its 
subtleties, since insertions of Ci renormalize the masses of the states in the effective theory 
and, therefore, modify the x-dependence of the "bare" states in (p^7|). In addition, there is a 
contact term arising from the action of the derivative on the ^-functions in the time-ordered 
product. We discuss these issues in Appendix C. One finds that the differences (Ei — E^) 
are in fact computable in terms of form factors introduced earlier. The relations are 



£ 7 - E' = A (v - v% A, + Vl (0o - Ai 



(4.29) 



E ya p — E' iaf3 



A(v- v') n A af3 - V^(v' a Jf3 - v'p^a) 4> 2 

+iv 1 a al3 (03 - A 2 • 



In particular, it follows that E[ = Ei for i = 3,6,9,10,11, and we will choose these five 
functions as a basis. Then a convenient way of writing the solution of the constraints 
imposed by the equation of motion is 



Ei+ w E 2 
Ei — E 2 
E[ - E' 2 



wE[ + E' 2 = E 3 , 
AAi +w(j)o, 
—AAi + 0o , 



E A + w E 5 
wE'^ + E's 

Ei + E 5 
E' 4 + E' 5 

E 7 + w E 8 
E 7 — E 8 
E' 7 - E' 8 



—Eq , 

Eq — E n 
1 

w + 1 
1 



E u + w<\) 2 - A(w - 1)A 2 



w 



1 



E n - 2 - A (w - l)A 2 



wE' 7 + E' t =E^ 

AA 3 + W 03 , 
-AA 3 + 03 , 



(4.30) 
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where we have introduced the nonsingular functions 



4> (w) - Ai^(w) 



w 



4> 3 (w) 



b 3 (w) - A 2 £(w) 



w 



(4.31) 



Consistency of the equations determining and E' 45 furthermore requires that, at zero 
recoil, 

2E 6 (1)-E U (1) = 2 (1). (4.32) 

The constraints imposed by the equation of motion allow us to prove a second theorem: 

Theorem 2: Matrix elements describing the mixed first order correc- 
tions to the current and to the Lagrangian vanish at zero recoil. 

It follows from the fact that, under the traces, 



E^(v,v) 



E X {1)+E 2 {1)-E 3 {1) 







E 



y a p(V,V)=l(T a i3 Vj 



EM) + E 8 (l) - E 9 (l] 







(4.33) 



with identical relations for E[. Thus, at zero recoil only genuine second order corrections 
to the current or to the Lagrangian contribute to hadronic form factors that are not kine- 
matically suppressed. The conservation of the vector current in the limit of equal masses 
then leads to relations between the universal functions which describe the corrections to 
the Lagrangian, and the parameters Ai and A 2 which, according to Theorem 1, describe the 
corrections to the current. These normalization conditions are the subject of the following 
subsection. 



D. Modified Wave Functions and Normalization 
Conditions at Zero Recoil 

In this section we have shown that at second order in the heavy quark expansion a total 
of 4 + 25 + 5 = 34 universal functions is necessary to parameterize, respectively, the effects 
of corrections to the current, of corrections to the effective Lagrangian, and of the combined 
corrections to both. The richness of the structures that arise might seem both impressive 
and frustrating, and the effort required to compute the various traces is quite considerable. 
However, only certain combinations of form factors appear in the final expression for any 
hadronic matrix element, and it is time to organize our results in a more transparent and 
convenient way by employing the concept of modified wave functions introduced in Sec. |TJ. 
The corrections proportional to 1/wq change the wave function for the initial state meson, 
but leave the final state unaffected (and vice versa for the terms proportional to l/nig,). 
Their effects can therefore be accounted for as in ( |3.12| ). On the other hand, the corrections 
proportional to l/rnqniQ' affect both mesons and can only be accounted for by a combined 
wave function. We can thus extend ( J3.16[ ) to second order by writing 
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(M'\ Q'TQ \M) 



-zli 2 zlif^w)tv{M'rM} 

' ■ tr{ M' T \P + (L* + -L + P (L M + -± 



2mr 



2m 

4m Q mQ 



2m c 

— tr{[(lf + JL?fW 
rig, I 1A + 2m Q / + / 4 

— tr{r 



If + 



2m Q 

1 -3M 



} 



2m, 



P+m™fPL + P-m MM 'P' 



-tI')P'_]TM} 

(4.34) 



MM 1 pi 



} + 



where we have performed the mass renormalization for the leading term. Here a "bar" de- 
notes Dirac conjugation combined with an exchange of velocities, polarizations, and masses. 
The virtue of ( |4.34| ) is that it allows an interpretation in terms of large and small compo- 
nents, reducing to a minimum the effort required to perform the traces. The structure of £ M 
is the same as that of L± in ( |3.13|) and involves six functions £i(w). The structure of m MM 
is more complicated and requires the introduction of twenty-four functions rrii(w). They are 
defined in Appendix A. 

Let us now discuss how the various second order corrections fit into this pattern. We 
start with the corrections to the Lagrangian, which according to ( [4.16[ ) and ( |4.17| ) preserve 
the P + projectors for the initial and final state. Hence, the fifteen universal functions P>i and 
Ci contribute to £^ only and appear in the three combinations £±,£2 and £3. Similarly, the 
functions Di contribute to m^ 1 and enter in the combinations m\ to m^. For the discussion 
of the corrections to the current we restrict ourselves to the operators in ( [4.1|) , which are 
obtained from tree level matching of QCD and HQET. As explained in Sec. OIL one can 
identify ip h with the small component of the full heavy quark spinor, and those terms lead 
to P projectors in the modified wave functions. The last operator in (|4. 1|) contains two such 
terms and consequently contributes to m^ 1 only. In fact, using the equation of motion its 
matrix elements can be written as 



(M'\ti (-ip)Tiph\M) = -tr{TP4^M4> al3 M'-? a }PL}. 



(4.35) 



By evaluating the bracket one readily computes the functions m§ to 77744, which appear in 
the parameterization of m__ . Because h'T^aVpG^h = —h'Tiv-Dip h, the other second 
order currents in (|4.1|) contain one small component and thus contribute to . To see this, 
we employ the equation of motion to write 

(M'| h'T la v p G af5 h\M) = -tr{7#rP..['f'Mv fi <f> af3 ]}. (4.36) 

The mixed corrections to the current and the Lagrangian have the same structure, since 

(M'\i fdxT{h'Ttph,Ci(x)} \M) = -tr{M'T P^ME^ + ^P+s^ME^p] }. (4.37) 



Thus, both fa and Ei enter in the functions £4, £5 and £q. Finally, the second matrix element 
in ( |4.24D determines the functions mis to m24, which appear in the decompositions of m^^f 
and m^ 1 '. 
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The complete set of expressions for l± and rrii is given in Appendix A, and in Appendix B 
we compute the meson form factors hi in terms of these functions. Let us now use these 
results to derive the normalization conditions which follow from the conservation of the 
vector current in the limit of equal heavy quark masses, tuq' = tuq. It implies that at zero 
recoil 



{M{v)\Q lfl Q\M{v)) = 2m M v fJ 



(4.38) 



for both pseudoscalar and vector mesons, which in terms of the meson form factors is equiv- 
alent to h + (l) = hx(l) — 1. It is now important that we have performed a mass renormal- 
ization in the first term in ( f4.34| ), since tum in ( f4.38| ) is the physical meson mass. Using the 
normalization of the Isgur-Wise function and Luke's theorem, we find from Appendix B in 
the equal mass limit 



Mi) 
Mi) 



1 + 
i 



im 2 Q 
1 



24(l)+mi(l)-m 8 (l) 



+ 



2£ 2 (1) + m 4 (l) + m 5 (l) - m n (l) - m 12 (l) 



+ 



Setting the coefficients of the second order terms to zero we obtain two conditions, which at 
tree level may be written as 



2B 1 (1) + 2d(l) + D 1 (l) - 3 2C 4 (1) + £> 4 (1) + 2C 5 (1) + D 5 (l) = -Ai 



(4.39) 



2fi 3 (l) + 2C , 3(l) + J D 3 (l)-2 



2C6(1) + D 5 (l)+2C 6 (1) + D 6 (i; 



A, 



More restrictive relations could be derived by including renormalization effects and requiring 
that the logarithmic dependence on txiq be the same on both sides of ( [4.391 ). The results of 
such an analysis will be presented elsewhere. 



V. APPLICATIONS AND SUMMARY 

Let us summarize the main results of our analysis. In total, thirty-four universal functions 
appear in second order of the heavy quark expansion of meson form factors. We have 
proved two theorems stating that, at zero recoil, the leading meson form factors do not 
receive contributions from mixed corrections to the current and the Lagrangian, and that 
the corrections to the current can be expressed in terms of Ai and A 2 . The number of 
universal functions is strongly reduced if one ignores radiative corrections and only considers 
the phenomenologically interesting cases of P — > P and P — > V transitions induced by a 
vector or axial current. Then all matrix elements can be parameterized in terms of l\ to £q 
and the five combinations (mi-mg), (m 2 + mg), (m^ — miQ), (mie + TOig) and (mn — mig). 
This can be seen from the relations given in Appendix B. 

In the following paragraphs we apply our results to semileptonic B decays and give 
estimates for some of the second order corrections. We also discuss the corrections to Luke's 
theorem, which arise at second order. For simplicity, we shall ignore radiative corrections. 
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A. Elastic Form Factors and B — > D iv Decays 



As pointed out in the introduction, the universal form factors of HQET describe the 
properties of the light degrees of freedom in the background of the color field of the heavy 
quark. From this point of view, the Isgur-Wise function is the elastic form factor that 
describes the overlap of the wave functions of the light degrees of freedom in the initial and 
final mesons moving at velocities v and v'. The normalization of £(w) at zero recoil reflects 
the complete overlap of the configurations of the light constituents in two infinitely heavy 
mesons with the same velocity. If finite-mass corrections are taken into account, the overlap 
decreases. In HQET the corresponding corrections are described by the functions L, and 
which represent the corrections to the wave function of a pseudoscalar (i — 1) or a vector 
meson (z = 2). At zero recoil, the first order corrections vanish, and using the expression 
for h+ and hi from Appendix B we obtain at second order 



(D(v)\ \B(v)) = 2^^ VfM {l + (e c - etfl^l) + ■■■}, 
(D*(v)\ \B*(v)) = -2y/m B *m D . e- e'* ^ {l + (e c - e b ) 2 £ 2 (l) + ■■■} 



(5.1) 



where Eq = l/2mg. 

In the nonrelativistic constituent quark model, the m^-dependence of the overlap integral 
comes from the mQ-dependence of the reduced mass of the light constituent quark, m^ ed = 
m q mQ/ '{rriQ + m q ). For an estimate of £i(l) we use the wave functions of the ISGW model 
[ 35| to obtain 



'1) = -3mi 



-0.75 GeV 2 . 



(5.2) 



For the numerical estimate we have identified the constituent mass of the light quark with 
A, since rriM = mq + m q in the ISGW model. 

The matrix element of the vector current between a B and a D meson enters the theoret- 
ical description of the decay rate for the semileptonic process B — > D iv. After contraction 
with the leptonic current a combination of the form factors h + and h_ appears [33], 



dw 



r<2 



cb\ 



487T 3 

x h + (w) 



m 3 D ( m B + "m D ) 2 (w 2 - 1) 3/2 



VSh_\ 



w ; 



(5.3) 



where S 



0.23 is the Voloshin-Shifman factor pj. At leading order in the 



heavy quark expansion the form factor is normalized at zero recoil, offering the possibility 
of a reliable determination of V c b for w ^ 1, provided that the corrections to the infinite 
quark mass limit are small. The first order power corrections are indeed suppressed by the 
Voloshin-Shifman factor and have been estimated to be ~ +2% [p5|j . Including the second 
order corrections, we find from Appendix B 



h + (l)-VSh-(V. 



l-(e c + e b ) S L 4 (l) + (e c - e 6 ) 2 [^(l) - 4(1)] + ±e c e b S A 



1 - 0.7% - 1.3% x 



A 2 



(5.4) 
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where we have used the heavy quark masses m c = 1.5 GeV and m& = 4.8 GeV, the con- 
stituent quark model estimate ( |5 . 2| ) , and the QCD sum rule results A « 0.5 GeV and 
-^4(1) ~ —A/3 [p5[| . For simplicity, the radiative corrections to h + and h- have been ne- 
glected. We conclude that, unless the coefficient £4(1) were unusually large, both the first 
and second order power corrections are small. Although not protected by Luke's theorem, 
the decay B — > D Iv thus allows for a reliable measurement of V c b. 
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B. Determination of V c f, from B — * D*£v Decays 

It has been observed in Refs. [0,^] that semileptonic B decays into D* vector mesons 
offer an almost model-independent measurement of V c b, since the 1/mg corrections to the 
decay rate vanish at zero recoil. In terms of the meson form factors one finds 

2 i \r. 12 



1 dr(E -f D*£v) G 2 F \V C 



lim - " — v ^ — = 1 V* 1 (m B - m^,*) 2 1^(1) T > ( 5 - 5 ) 



V^ 2 — 1 dw 47T 3 



and /iAi(l) is protected by Luke's theorem [15). Thus the determination of from an 



extrapolation of the spectrum to w ^ 1 is model-independent up to terms of order 1/m 2 . 
From Appendix B we obtain at second order 



where 



h Al (l) = 1 + (e c - e 6 ) [e c € 2 (l) - Mi(l)] + £^A, (5.6) 



A = £ 1 (l)+£ 2 (l)+m 2 (l)+m 9 (l) 

Ai + 2 A 2 + 4 Ld 4 (1) + 2D 5 (1) + D 6 (l)l • (5.7) 



_ 4 
3 

Using ( |5.2| ), the first correction in ( p.6| ) is estimated to be — 3A 2 (£ C — £ b ) 2 w —3.9%. Con- 
cerning the second term we observe that, except for Ai and A 2 , the coefficient A depends 
only on form factors which arise from a double insertion of the chromo-magnetic moment 
operator of C\ in ( |2.4|) . We shall argue below that these terms are expected to be very small. 
Neglecting them, and using fl2.15f) as well as the sum rule estimate Ai ~ 1 GeV 2 [[17 
obtain e c e h A w 5.7%, and thus 



wc 



/i jtl (l)-l«2%. (5.8) 

The main uncertainty in this estimate arises from the uncertainty in Ai, as discussed at the 
end of Sec. |T[ In the extreme case Ai = we would obtain Ka-SX) ~ 1 ~ — 3% instead 



of (5.8). However, in any case the second order correction is small because of a partial 
cancellation of the two terms in (|5.6|) , suggesting that the theoretical uncertainty in this 
method of extracting V c b is less than a few percent. 



It has been claimed in Ref. [36] that QCD sum rules would predict a second order 
correction to ft^ii(l) of as much as — 10%.! In view of our estimate ( |5.8| ) this assertion seems 
unacceptable. Even for Ai = it would imply that ^i(l) and ^ 2 (1) would have to exceed the 
quark model prediction (|5.2|) by a factor of three. 



C. Second Order Corrections to Luke's Theorem 

At the end of Sec. |TTT], we discussed the fact that Luke's theorem protects the meson 
form factors h + , h Al , h\, and h 7 from first order power corrections at zero recoil. Although 



2 It has been pointed out in Ref. [25| that the argument given in Ref. [Q has no theoretical 
foundation. 



22 



our results show that there is no such nonrenormalization theorem at second order, the 
structure of the 1 / m 2 corrections to these four form factors is particularly simple and allows 
for a semi- quantitative estimate. At zero recoil, the expression for hr is 

h 7 (l) = 1 + (e c - e b ) 2 £ 2 {\) + e c e b A' (5.9) 

with 

A' = | Ai - 2 A 2 + 4 [D 4 (l) - As(l)] • (5.10) 

The other three form factors have been given in ( |5.1Q and ( |5.6|) . We observe that there 
is always a correction involving £i(l) or £2(1), depending on whether one deals with a 
pseudoscalar or a vector meson, respectively. Using the quark model estimate ( |5.2j ), this 
term becomes approximately —4%. It smallness naturally results from the squared difference 
(e c — Eb) 2 - In addition, for and I17 there is a term proportional to £ c e&, which depends 
on the mass parameters Ai and A2 as well as on form factors arising from a double insertion 
of the chromo- magnetic moment operator. Neglecting these latter terms, this correction can 
be estimated based on a model calculation of Ai, since A2 is known from the experimentally 
observed mass splitting between vector and pseudoscalar mesons. QCD sum rules predict 
that Ai is positive, and the corresponding correction tends to cancel the terms proportional 
to £i, which are negative. As a result, the form factors and h-j can only receive small 
1/m 2 corrections at zero recoil. 



D. Limit of Vanishing Chromo- Magnetic Interaction 

Detailed QCD sum rule analyses of the universal functions that appear at order 1/m 
in the heavy quark expansion show that the form factors A2 and A3, which arise from the 
insertion of the chromo- magnetic moment operator in C\, are much smaller than the other 



two functions, A\ and £3 [25.37]. The coarse pattern of the 1/m corrections can be well 
described by setting A2 and A3 to zero, corresponding to the fictitious limit of vanishing 
field strength, G a ^ — > 0. Let us see what kind of simplifications the same approximation 
implies at order 1/m 2 . 

We start with the corrections to the current. In the limit G a ^ — > the functions <pi, 
4>2 and 03 vanish, and according to ( |4. LQ ) this implies the vanishing of Ai and A2. It then 
follows that 0o = (w — 1) <f) vanishes at zero recoil, and all corrections to the current can be 
described by the single function 0. Similar simplifications occur for the corrections to the 
Lagrangian. Here all universal functions except Bi,C±, and D\ vanish in the limit G af3 — > 0. 
At tree level, one obtains from ([4.39|) the zero recoil condition 



Bi(l) + d(l) + \ D X {1) = . (G a ? - 0) (5.11) 

Finally, the combined corrections to the current and to the Lagrangian are entirely param- 
eterized by the form factor E3, since E 6 , Eg, E w , and £ u vanish in the limit of vanishing 
field strength. 

In the fictitious limit of vanishing chromo-magnetic interaction, the set of thirty-four 
universal form factors is thus reduced to only five functions, a combination of which vanishes 
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at zero recoil. Although we are aware of the fact that such an approximation can only give 
us a very simplified picture, we still believe that it might be useful for an analysis of the 
structure of the dominant terms. The expressions arising for the functions £i and m; in this 
limit can readily be obtained from the general formulas given in Appendix A. 



E. Summary 

Using the heavy quark effective theory, we have performed the expansion of matrix ele- 
ments of heavy quark currents between pseudoscalar or vector mesons up to second order in 
inverse powers of the heavy quark masses. The general description of the power corrections 
arising at order l/m 2 involves a set of thirty- four Isgur-Wise form factors, which are univer- 
sal, mQ-independent functions of the kinematic variable w = v ■ v ' . These form factors are 
defined in terms of matrix elements of higher dimension operators in the effective theory. 

Apart from some normalization conditions imposed by vector current conservation, the 
universal functions are hadronic quantities which cannot yet be predicted from first prin- 
ciples. Nevertheless, we have argued that in certain cases of phenomenological interest 
the l/m 2 corrections are parameterically suppressed. In particular, the corrections to the 
semileptonic decay rates for B — > D £ v and B — > D*£ v at zero recoil are estimated to be 
small, not exceeding a few percent. Our results thus support the usefulness of the heavy 
quark symmetries for an accurate determination of the weak mixing parameter V c b from 
these decay modes. 

Although the structure of second order corrections to various decay rates is quite complex, 
we believe that a classification in terms of universal form factors is still a useful concept. In 
particular, this might provide a framework in which to analyze various models. For instance, 
we have shown that the second order corrections to elastic form factors arising from the ttiq- 
dependence of the reduced mass of the light constituent quark in a nonrelativistic quark 
model are accounted for by our functions l\ and £ 2 , and an estimate of the effect gives 
£i(l) m ^2(1) ~ — 0.75GeV 2 . This information can then be used to predict corrections to 
other form factors, whose dependence on £\ and £2 is known from heavy quark symmetry. We 
have also suggested that, for an estimate of the dominant corrections, one might consider 
the limit of vanishing chromo-magnetic interaction, in which only five of the thirty-four 
universal form factors remain. The usefulness of such an approximation is supported by 
QCD sum rule calculations of the form factors appearing at order l/m in the heavy quark 
expansion. 

The analysis presented here for mesons can straightforwardly be extended to other 
hadrons containing a single heavy quark. The particularly interesting case of the A baryons 



is discussed in Ref. 27 
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APPENDIX A: COMPUTATION OF THE MODIFIED WAVE FUNCTIONS 



According to ( |3.13| ), the general structure of the modified wave functions £± introduced 



in flOlj) is 



op 

-+ 
9 p 

iV 



m(- 7b )*i, 

\f {£ 2 + e-v'£ 3 

M(- 75 )4, 

£ 5 + e-v'e 6 



(Al) 



The coefficients £\ are functions of w = v -v ' . Similarly, we choose the following decomposition 
for the product wave functions m+J^ 1 : 



m 
m 
m\ 
rn 



pp 
++ 

PV 

++ 

.VP 

++ 

vv 

++ 



VMM' mi (-75) 75 = -VMM' mi , 
VMM' [m 2 (-75) f* + m 3 (-75) e'*-v 

Vmw 1 -' j - 1 ' 1 -- py 



mo ^75 + 777.3 6-1/75 



++ ' 



"2 

VMM' m 4 +m 5 e- e* + m 6 / e'*- ?; + m' 6 ^* e • 7/ + m 7 e ■ v' e 1 *- v 



(A2) 



An identical decomposition with functions ms to 77744 applies for m__ . Finally, we define 
,pp 

= - V 1V11V1- 7/615 , 

_i (A3) 



rn. 



777 



PV 



-VMM' 77715 , 

VMM' [m 16 (-75) <* + m 17 (-75) e'*^ 



7/7 



VP 



VMM' 



nii % (^ + mi 9 e-v 75 



.vv 



VMM' f 77720 / + 77721 6 • e'* + 777 22 / €*■ V + 777 23 € • l/ + 777 24 € • l/ e'*- V 



and m^f^ is described by a set of related functions m' { , since 777 ^ f H M ' = m+'J* 1 . In these 
expressions a "bar" means Dirac conjugation combined with an exchange of velocities, po- 
larizations, and meson masses. Also, because of radiative corrections the functions £{ and 
77ij depend logarithmically on the heavy quark masses, and 777^ are related to by an 
interchange of mq and 777q/ in the renormalization factors. At tree level there is no such 
difference, and £{ and rrii are universal, mg-independent functions. 

The tree level expressions for these functions can be obtained by evaluating the various 
traces, as explained in Sec. fV[ We find: 

£1 = (Ai + 3A 2 ) £ + Bi + 2(w - 1)B 2 + 3v3 3 + C x + 2{w - 1)C 2 + 3C 3 
-3C 4 - 9C 5 - 6C 6 + 2{w 2 - 1)(2C 7 + C 8 ) - 4(w - 1)(C 9 + C 12 ) 

£2 = (Ai - A 2 ) f + Bi - B 3 + Ci - C 3 - 3C 4 -C 5 + 2C 6 
+2(w - 1) [(w + 1)C 8 -C 9 - C 10 + 3C U + C 12 

£ 3 = -2B 2 - 2C 2 + 4(w + 1)C 7 + 2C 9 - 2Ci - 10C n - 10C 12 

£ 4 = -ALi - w(0 o + 30 3 ) -(w+ l)0i + 402 + 30 3 
-2E 3 - A(w + 1)E 6 - 6E g + 2(w + l)E w - 4E n 
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4 



w{(p 



(w + l)0i + 20 2 + 3 + 2(w + 1)E 10 - 2E 



li 



w+l 



AL 2 + ±(w - 1)L 3 + w(0 o - 3 ) + + 1)01 - 02 
- £ 3 + 2(w + l)E Q + E 9 - 2{w + 1)E 10 + E u 



m 1 = D 1 + 2(w - l)D 2 + 3L> 3 - (2w + 1)D 4 - 9L> 5 - 6L> 6 

-2(w - 1) [{w + 1)(2D 7 + D 8 ) - 4D 9 - 8D10 
m 2 = D x + (w - l)D 2 + D 3 + D 4 + 3D 5 + 2D 6 - 2(w - l)(D 9 + D w ) 
m 3 = -2D 2 + 2D A + 2(w + 1)(2D 7 + D 8 ) - 2D 9 - 10D 10 , 
m 4 = D x - D 3 + (2w - 1)D 4 - D 5 - 2{2w - 1)D 6 

+2(w - 1) [(w + 1)D 8 - 2D 9 + 2D 10 " 

m 5 = -4wD A + 4{2w - 1)D 6 - 4{w - 1) [{w + 1)D 8 - 2D 9 + 2D l0 
m 6 = -2D 2 - 2D A + 4D 6 - 2{w + 1)D 8 + 6D 9 - 2D 10 
m 7 = 4L> 4 - 8D 6 - 4D 7 + 4wD 8 - 8D 9 + 8D 10 



i 

w+l 



1 

w+l 



302 + |03 1" lJ - 



m 8 
m 9 

77lio 

mn 

"7l2 

™i3 = -0 + ^tt 02 + i03 + A(L 4 - 2L 5 ) 

77li4 



0o + ^103 + ^ [(™ + l)0i - 602 - 2AL 4 

w + 1)01 + 2(2 - u>)02 + 303 - A(w - 1)L 4 
AL 



-00 - (W + 1)01 + 202 + 203 

20o + 2u>0i - [2w0 2 + (2m + 1)0 3 + A(w - 1)L 5 



w+l 

' w r0 + 



w+l 



(w+l) 2 



4{w + 1)01 - 2(3u> + 4)02 + {w + 2)03 
+ 4AL 4 - 2A(4 - w)L 5 



m 15 

mi 6 
777,17 

77718 
77719 

m 20 
m 21 
m 22 

m 23 

m 2 4 



4E 



ii 



— ALi + 0o + 303 - 202 
-2E 3 - 4{w + 1)E 6 - 6E 9 + 2(w + l)E w 
-AL 2 + 0o - 03 - 2E 3 + 2E 9 + 2E n 

l)L 3 - 02 - 2{w + 1)E 6 + (w + 1)E 10 + En 

00 + 303 - 202 



w+l 

-ALi 



ALi - O 



,^ 30 3 + 2 2 - E 3 - 2(w + 1)E 6 - 3E 9 

-AL 2 + 0o - 03 - 2{w + 1)E 10 + 2E U 
4(w + 1)E 10 - 4E n 

A(w - 1)L 3 - 202 - 2(w + 1)E 10 + 2E n 
- AL 2 + 0o - 03 + E 3 - E 9 - (w + 1)E W 



(w 



w+l 



2 

w+l 



U 2 -! 



;(w - l)L 3 - 02 + (w + 1)E 6 + (w 2 - l)Eio - wE u 



l)E 



10 
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Eq. ([4.32 ) ensures that there is no pole in m 24 as w — ► 1. Note that the first term in 4 and 
4 compensates the mass renormalization performed in ( |4.34| ). 



APPENDIX B: MESON FORM FACTORS 

Let us set Eq = l/2rriQ. Then to second order in the heavy quark expansion the meson 
form factors hi introduced in (|3.17 ) are given by: 



h + = i + (e c + e b ) Li + (e 2 c + e 2 b ) t x + e c e b {nix - m 8 ) 



h. 



(e c - e b ) L 4 + (e 2 c - e 2 b ) 4 



(Bl) 



h v = £ + e c (L 2 - L 5 ) + e b (L x - L 4 ) 

(£ 2 - 4) + £& (4 - 4) + ^c^b [(™<2 + mg) - (miB + mis)] 



£ + e c (L 2 - ^q-[L 5 ) + £ 6 (Li - {j^L, 



+el (4 -^4) + ^4 -^4 )+■--,-.-/, 



^m 2 + m 9 ) - ^pj(mi 6 + mi 8 ) 



e c (L 3 + L 6 ) + e 2 c (4 + 4) + e c £: 6 [(m 3 - m 10 ) - (m 17 - m 19 )] 
£ + e c (L 2 - L 3 - L 5 + L 6 ) + £ 6 (Li - L 4 ) 

+e 2 c (4-4 -4 + 4) + ^(4 -4) 

+e c e 6 [(m 2 + m 9 ) - (m 3 - m 10 ) - (m 16 + m 18 ) - (m 17 - mi 9 )] 



(B2) 



(B3) 



hi = £ + (e c + e 6 ) L 2 + (ej! + e 2 b ) £ 2 + e c £: 6 [(m 4 - m n ) + (m 5 - m 12 )] 
^2 = (e c - e&) ^5 + (£c - e?) 4 

/i 3 = £ + £c [L 2 + (w - 1)L 3 + L 5 - (w + 1)L 6 ] + e b (La - L 5 ) 

+^ [4 + ( w - 1)4 + 4 _ („, + 1)4] + el (4 _ 4) 

+e c £: 6 [(m 4 - mn) + (m 5 - m i2 ) - (w - l)(m 6 - mi 3 ) - (w + l)(m 22 + m 23 )] 
^4 = h 3 (e c <-> e 6 ) 

/i 5 = e c (L 3 - L 6 ) + (4 - 4) + £ c £b [{m 6 - m 13 ) + (m 7 - m i4 ) - (m 22 + m 23 )] 

= h(e c «-> £ 6 ) (B4) 



^7 = f + (e c + ^2 + (^c + £ b) 4 + £ c £fe (m 4 - m u ) 
h 8 = h 2 



(B5) 



These relations are valid at tree level. The radiative corrections to the leading and subleading 
terms in the heavy quark expansion have been calculated in Refs. |3|,[38|,[3!| . 



APPENDIX C: MODIFIED WARD IDENTITIES 

Here we derive Ward identities which relate the derivative of the matrix elements in 
( p.lO| ) to the matrix elements in ( |4.24| ), in which a derivative acts on the current. These 
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identities are needed in Sec. |TV|C to express the universal functions E[ in terms of and 
other form factors. Let us consider the following matrix element: 



(M'\J(z)\M + 6M) = (M'\J(z)\M) + (M'\idxT{J(z),Cx(x)}\M). (CI) 

2m Q J 

J(z) is a heavy quark current in the effective theory, \M) is an eigenstate of £hqet, and 
\M + 6M) denotes an eigenstate of £hqet + ^jA- ^ n con t ra st to (p.7|) , we have 



\M + SM) Z = exp 



-i A + 



Am 2 M 
2mn 



v ■ z 



\M + 6M) . 



(C2) 



Using this fact, we find to order 1/tuq 

id* (M'\ J(z) \M + 6M) = id* (M'\ J{z) \M) 



Am 



M 



2m, 



w 7 (M'\ J(z) \M) 



A 



2m f 



(v - v% (M'\ i dxT{ J(z), Ci(x) } \M) . 



(C3) 



Collecting terms of order 1/mq, we thus obtain 



id* - A (v - v% (M'|i/da;T{ J(z),£i(x) } |M) = Am> 7 (M'| |M) . (C4) 



On the other hand, carrying out the derivative acting on the time-ordered product gives 



id* (M'\i dxT{ J(z), d(x) } \M) = (M'\ i / dxT{ id y J(z), d(x) } \M) 



(C5) 

-v 7 (M'| /i' T P + [{iDf + Z s a/3 G Q/3 ] h \M) . 
Combining ( |C4[ ) and ( |C5| ) , we find for the universal form factors the relations given in ( [4.291 ). 
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FIGURES 



FIG. 1. Diagrams representing the first order power corrections to meson form factors in HQET: 
(a) corrections to the current, and (b) corrections to the effective Lagrangian. The squares represent 
operators of order 1/mg or 1/rriQi. 

FIG. 2. Diagrams representing the second order power corrections to meson form factors in 
HQET: (a) corrections to the current, (b) corrections to the effective Lagrangian, and (c) mixed 
corrections to the current and the effective Lagrangian. The black squares represent operators of 
order 1/rriQ or l/mgi, the open ones denote operators of order 1/rriQ, 1/rriQ,, or 1/mQrriQi. 
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